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1.  INTRODUCTION 


1 . 1  Background 

Nations  which  have  signed  the  Partial  Test  Ban  Treaty  have  agreed  not  to 
conduct  nuclear  weapons  tests  in  the  atmosphere,  underwater  or  in  space. 
Nations  which  have  signed  the  Nuclear  Non-proliferation  Treaty  have  agreed  not 
to  develop  nuclear  weapons  or  assist  others  to  do  so. 

Monitoring  the  environment  to  detect  nuclear  weapons  tests  is  of  interest 
for  several  reasons.  It  is  desirable  to  assure  that  there  is  no  violation  of 
the  Partial  Test  Ban  Treaty  or  Non-proliferation  Treaty  by  signature  nations, 
and  that  non-signature  nations  have  not  developed  nuclear  devices. 

The  technical  data  on  which  monitoring  systems  are  based  has  been 
obtained  over  a  research  period  of  four  decades.  The  data  has  been  obtained 
from  nuclear  tests  performed  prior  to  the  inception  of  the  various  applicable 
treaties,  from  legal  nuclear  tests,  from  nuclear  tests  performed  in  banned 
environments  by  nations  which  have  not  signed  the  applicable  treaties,  from 
simulation  tests  performed  with  chemical  high  explosives,  and  from  theoretical 
investigations  of  the  nuclear  explosion  processes. 

The  technical  basis  for  detection  of  nuclear  detonations  in  the  banned 
environments  using  remote  sensors  is  well  established.  In  addition  to  the 
generation  of  radioactive  debris,  a  nuclear  detonation  in  the  atmosphere  is 
characterized  by  a  double  light  flash  and  radiation  which  affects  the 
ionosphere.  For  an  underwater  detonation,  high  level  acoustic  signals  are 
generated  and  the  efficiency  of  seismic  wave  production  is  very  high.  The 
technical  problems  of  monitoring,  if  any,  stem  from  implementation  to  obtain 
sufficient  information  in  a  timely  fashion. 
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It  has  been  suggested  that  we  can  be  certain  that  a  nuclear  event  has 
occurred  only  if  radioactive  debris  is  found  at  the  location  predicted  on  the 
basis  of  other  observations.  This  is  a  rather  rigid  requirement  which  may  be 
difficult  to  satisfy.  If  the  event  occurs  in  a  remote  oceanic  region  under 
circumstances  where  there  is  no  radioactive  debris  deposited  on  land,  the  time 
delay  between  the  initial  observation  and  the  arrival  at  the  scene  of  vehicles 
equipped  with  proper  scientific  equipment  may  preclude  the  acquisition  of 
debris  prior  to  its  natural  dispersal. 

Without  comment  as  to  its  wisdom  or  reasonableness,  we  adopt  the  view 
that  a  viable  monitoring  system  must  provide  a  capability  for  satisfying  the 
requirement  for  retrieval  of  radioactive  debris.  The  total  system  capablity 
can  then  be  broken  down  into  three  broad  categories: 

a.  A  remote  sensor  system  capable  of  monitoring  the  known  environmental 
effects  with  real-time  detection  and  localization  of  suspect  events. 

b.  A  signal  processing  system  capable  of  analyzing  the  data  obtained  for 
the  suspect  event,  confirming  the  detection,  improving  the 
localization,  and  estimating  the  yield  and  other  pertinent  detonation 
parameters. 

c.  A  capability  for  rapid  deployment  of  an  appropriate  vehicle  to  the 
detonation  site  for  retrieval  of  radioactive  debris. 

In  this  study,  we  will  assume  that  the  first  of  these  system  categories  - 
the  sensor  system  -  is  operable.  The  extent  to  which  appropriate  vehicles  for 
retrieval  of  radioac  ive  debris  are  maintained  is  part  of  a  more  general 
overall  readiness  problem  and  will  not  be  addressed  here.  We  will  consider 
the  second  category  restricted  to  underwater  or  over-ocean  atmospheric 
detonations  and  the  requirements  this  may  impose  on  the  sensor  systems. 


1 . 2  Hydroacoustic  Sensor  System 

We  assume  that  a  total  monitoring  system  has  already  been  deployed  which 
includes  the  following  components  and/or  capabilities: 

a.  Hydrophones  are  located  in  each  ocean  basin  of  interest  to  detect  the 
hydroacoustic  signals  generated  by  underwater  or  above-water 
detonations,  and  possible  T-phase-1 ike  signals  generated  by 
underground  detonati  ns  located  near  a  coast. 

b.  A  seismometer  network  for  detecting  the  seismic  signals  generatea  by 
underground  or  underwater  detonations. 

c.  A  satellite  system  for  detecting  the  double  light  blip  generated  by 
an  atmospheric  detonation. 

It  is  further  assumed  that  appropriate  elements  of  the  system  are 
monitored  in  real  time  to  detect  and  identify  suspect  events.  When  a  suspect 
event  is  identified,  the  hydrophone  output  recordings  would  be  carefully 
examined.  If  direct  arrivals  are  received  at  three  hydrophone  locations,  the 
geographical  position  of  the  detonation  can  be  determined.  If  direct  arrivals 
are  not  observed  on  three  hyonphones,  other  supporting  information  is 
necessary.  For  an  atmospheric  detonation,  this  could  be  the  detonation  time 
based  on  the  satellite  observation  of  the  time  of  occurence  of  the  double 
light  blip.  In  any  event,  we  assume  that  procedures  for  determining  the 
geographical  location  of  the  detonation  based  on  tri angulat ion  using  direct 
arrival  times  and  supplementary  information  have  already  been  implemented. 

It  is  perhaps  more  realistic  to  recognize  that  direct  signals  will  rarely 
be  received  on  three  hydrophones.  This  is  due  to  a  combination  of  factors. 
First,  the  ocean  area  which  has  to  be  monitored  is  very  large  while  the  number 
of  hydrophones  employed  is  necessarily  limited.  Second,  a  nation  which  wishes 
to  evade  detection  in  performing  a  test  with  an  underwater  or  above-water 
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detonation  can  take  advantage  of  topjgraphic  shielding  to  greatly  reduce  the 
probability  that  a  direct  signal  will  be  received,  even  if  it  does  not  know 
where  the  monitoring  hydrophones  are  located.  However,  it  is  much  more 
difficult  to  select  a  detonation  location  for  which  basin  reverberation  is  not 
received.  It  has  been  found  that  the  reverberation  received  on  a  single 
hydrophone  is,  in  principle,  sufficient  to  permit  the  application  of 
triangulation  principles  to  determine  the  geographical  location  of  the 
source.  The  techniques  for  doing  this  are  in  the  very  early  stages  of 
development.  Nevertheless,  the  use  of  reverberation  as  a  means  for 
determining  the  location  of  the  detonation,  and  estimating  the  yield  and 
depth/ al  titude ,  is  a  procedure  which  should  be  considered  as  an  addition  to 
the  total  system  capability.  However,  the  major  thrust  of  this  study  is  the 
definition  of  procedures  for  determining  the  yield  and  depth/altitude  of  the 
detonation  from  an  analysis  of  the  direct  arrival  observed  on  a  hydrophone 
once  the  geographical  location  of  the  detonation  is  known,  without  our 
specifying  how  the  location  is  determined. 


Chapter  2  is  divided  into  three  major  sections.  The  first  section 
presents  the  computational  scheme  for  determination  of  yield  vs  depth/altitude 
for  a  suspect  underwater  detonation  given  that  an  underwater  acoustic  signal 
has  been  received  and  appropriately  processed,  and  that  the  location  of  the 
suspect  event  is  known.  The  second  section  is  a  discussion  of  the  method  to 
be  employed  for  processing  of  the  received  underwater  acoustic  signal.  The 
final  section  identifies  potential  sources  of  error  in  the  determination  of 
yield  vs  depth/ altitude  and  provides  a  brief  quantitative  analysis  based  on 
assumed  error  statistics. 


Chapter  3  provides  a  brief  summary  of  the  procedures  and  discussions 
presented  in  Chapter  2. 

Chapter  4  presents  the  authors'  recommendations  resulting  from  the  Phase 
1  effort. 

Appendix  A  provides  a  discussion  of  the  methodology  to  be  employed  for 
determination  of  the  geographic  location  of  a  detonation  based  on  the 
pressure-time  history  of  the  shot-generated  reverberation  received  at  a  single 
hydrophone.  Also  included  are  brief  descriptions  of  two  extant  computer- 
implementea  reverberation  time-history  models. 

Appendix  B  presents  a  brief  discussion  of  the  detectability  of  seismic 
signals  generated  by  over-water  nuclear  detonations.  A  method  for  relating 
the  seismic  energy  generated  by  deep  underwater  and  atmospheric  detonations  is 


2.  COMPUTATIONAL  SCHEME 


Conceptually,  the  software  system  for  the  determination  of  the  yield  and 
depth/altitude  of  a  nuclear  event  can  be  divided  into  two  parts.  One  part 
comprises  those  algorithms  which  are  employed  to  determine  a  yield-vs- 
depth/ al  titude  curve  for  such  an  occurrence  based  on  the  received 
hydroacoustic  signal.  The  determination  of  this  functional  relationship  is 
the  heart  of  the  classification  (yield  and  depth/altitude)  process  and  is 
given  major  emphasis  in  this  study.  Its  determination  involves  the 
application  of  nuclear  blast  theory  and  basic  acoustic  principles.  The  other 
part  of  the  system  comprises  those  algorithms  employed  to  determine  where  the 
event  lies  on  the  previously  determined  yield-vs-depth/al titude  curve.  This 
part  of  the  process  brings  to  bear  the  capabilities  of  the  total  monitoring 
system  in  providing  clues  as  to  the  nature  of  an  explosive  event. 
Determinations  based  on  the  information  provided  by  the  overall  system  will  be 
quasi-analytical ,  will  necessarily  involve  a  certain  degree  of  subjectivity 
and,  in  all  likelihood,  will  never  be  sufficient  to  pinpoint  an  event  on  the 
yield-vs-depth/al titude  curve.  The  best  that  can  be  hoped  for  are 
specifications  of  yield  and  depth/ altitude  within  such  bounds  as  to  be 
significantly  useful  to  those  at  the  decision  making  levels  of  government. 

2.1  Determination  of  Yield  Vs  Depth/Altitude 

The  following  discussion  is  prefaced  by  these  basic  assumptions.  First, 
it  is  assumed  that  the  geographic  location  of  the  detonation  has  been 
determined.  Secondly,  the  hydrophone  is  assumed  to  be  located  within  the  deep 
sound  channel  in  water  of  sufficient  depth  to  support  deep  RSR  propagation. 
Thirdly,  it  is  assumed  that  there  is  no  intervening  land  mass  along  the  great 


circle  path  between  the  detonation  site  and  the  hydroacoustic  sensor, 
precluding  reception  of  a  direct  signal.  And,  finally,  the  distance  between 
the  detonation  site  and  the  receiver  is  assumed  to  be  sufficiently  great  that 
the  geometric  spreading  of  the  acoustic  energy  at  the  separation  distance  is 
essentially  cylindrical.  It  is  important  to  note  that  the  system  places  no 
restriction  on  the  propagation  conditions  at  the  detonation  site.  That  is, 
the  water  column  in  the  vicinity  of  the  site  may  or  may  not  support  RRR  or  RSR 
propagation. 

Theoretical  considerations  based  on  nuclear  detonation  phenomenology 

dictate  that,  from  a  computational  standpoint,  we  divide  the  depth/ al t itude 

axis  into  three  regions  -  underwater,  atmospheric  and  near-surface.  For  a 

given  detonation  yield,  the  underwater  region  refers  to  depths  greater  than 
1/3 

21 W  feet,  where  W  is  the  yield  in  kilotons.  For  detonation  depths 

1/3 

greater  than  21W  feet,  direct  application  of  usual  hydroacoustic  methods 

1/3 

may  be  employed.  For  underwater  source  depths  shallower  than  21W  feet, 
all  of  the  water  above  the  detonation  is  vaporized  and  new  phenomena  must  be 
considered.  The  theory  developed  in  Reference  1  relevant  to  the  determination 
of  a  yiel d/al t itude  curve  for  atmospheric  shots  is  predicated  on  the 
assumption  that  the  detonation  occurs  at  such  height  that  the  shock  wave 
separates  from  the  fireball  (hydrodynamic  separation)  before  interacting  with 
the  sea  surface.  In  the  atmosphere,  hydrodynamic  separation  occurs  at  a  range 
of  15W^/3  feet,  where  W  is  the  yield  in  kilotons.  The  atmospheric  region, 
then,  is  that  portion  of  the  axis  above  lSW1^  feet.  Thus,  the  near-surface 
region  extends  from  21W^^  feet  in  water  to  15W^/3  feet  in  air. 
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2.1.1  Yield  Vs  Depth  -  Underwater  Case 

The  first  step  in  the  specification  of  the  yield-vs-depth  curve  for  an 
underwater  shot  is  the  determination  of  a  source  level-vs-depth  curve  based  on 
the  received  acoustic  signal.  The  determination  of  source  level  vs  depth 
requires  that,  for  the  analysis  frequency,  the  detonation  site-to-receiver 
propagation  loss  be  estimated  for  detonation  depths  throughout  the  water 
column.  The  source  spectrum  level  at  each  depth  is  given  by 

SL  =  RL  +  PL  (1) 

where 

SL  is  the  source  spectrum  level 

RL  is  the  received  signal  spectrum  level 

PL  is  the  predicted  propagation  loss 

Thus,  the  source  level  determination  at  any  given  assumed  source  depth  is  that 
level  which  would  produce  the  actual  received  signal  level  if  the  detonation 
had  occurred  at  the  given  depth.  Because  the  classification  system  is 
intended  for  use  in  situations  where  the  detonation  site-receiver  separations 
are  quite  large,  and,  therefore,  where  uncertainties  in  the  detail  of  the 
propagation  loss-vs-range  function  are  large,  broad  range  averaging  -  say 
30  nm  -  will  be  employed  in  the  propagation  prediction  process.  The 
specification  of  an  appropriate  transmission  loss  model  (or  models)  for 
incorporation  into  the  classification  system  is  beyond  the  scope  of  this 
study.  Suffice  it  to  state  here  that  the  prediction  requirements  can 
certainly  be  met  through  the  use  of  one  or  more  of  the  computer-implemented 
models  currently  available  to  the  Navy. 

The  next  computational  step  is  the  conversion  of  the  source  level-vs- 
depth  curve,  as  determined  above,  to  the  corresponding  yield-vs-depth  curve. 

If  the  actual  burst  occurs  in  the  atmosphere  or  at  such  depth  in  the  ocean  as 


to  result  in  the  venting  of  the  bubble  formed  by  the  detonation,  the  received 
underwater  signal  is  the  result  of  acoustic  energy  generated  solely  by  the 
shock  wave.  If  the  burst  occurs  underwater  below  venting  depth,  the  received 
signal  reflects  the  sum  of  the  energy  generated  by  the  shock  wave  and  bubble 
pulse.  However,  the  bubble  pulse  period  for  nuclear  charges  is  so  long  that 
the  energy  contributed  by  the  bubble  is  limited  to  very  low  frequencies.  If 
we  restrict  the  analysis  to  frequencies  greater  than  about  10  Hz,  the  bubble 
energy  is  negligible  and  the  received  signal  can  be  assumed  to  have  been 
generated  solely  by  the  shock  energy.  Yield  determinations  can  therefore  be 
made  using  acoustic  source  level-vs-yiel d  data  based  on  the  nuclear  shock  wave 
only.  In  this  document,  an  analysis  frequency  of  20  Hz  is  assumed. 

For  detonations  underwater,  it  can  be  shown  that  nuclear  shock  wave 
source  levels  vary  approximately  as  7.5  log  W.  Expressed  another  way,  source 
levels  change  by  7.5  dB  per  order  of  magnitude  change  in  yield.  This 
relationship  permits  one  to  determine  -  for  a  given  frequency  -  the  yield  of  a 
detonation  based  on  its  source  level  and  a  known  yield  and  source  level  used 
as  a  reference.  The  algorithm  is  given  by 


where 

W  is  the  yield  to  be  determined 

WREF  the  reference  yield 
SL  is  the  source  level 

^■REF  the  reference  yield  source  level 
Based  on  information  provided  in  Reference  2,  it  is  found  that  for  a 
1-kiloton  detonation  the  source  level  at  20  Hz  is  268.4  dB//gPa/Hz  at 
1  yard.  Using  1  kiloton  as  the  reference  yield,  Eq.  (2)  becomes 


W  *  10 


(3) 


where 

W  is  the  yield  in  kilotons 
SL  is  the  source  level  in  dB//yPa/Hz  at  1  yard 
Eq.  (3)  is  the  algorithm  to  be  used  in  determining  the  yield-vs-depth  curve 
for  the  underwater  case.  The  curve  is  terminated  at  that  depth  above  which, 
for  any  given  source  depth,  the  value  of  21W^  is  greater  than  the  given 
depth. 


2.1.2  Yield  Vs  Altitude  -  Atmospheric  Case 

For  an  over-water  nuclear  detonation,  there  are  two  wave  types  which 
can  penetrate  the  air-water  interface  and  propagate  to  long  range  via  RSR 
paths.  These  are  the  geometric  wave  and  the  lateral  wave.  The  geometric 
wave  is  that  wave  whose  propagation  can  be  described  by  normal  ray  acoustic 
methods.  The  waterborne  lateral  wave  results  from  the  fact  that 
exponentially  decaying  inhomogeneous  plane  waves  present  in  the  expansion  of 
the  spherical  wave  issuing  from  the  detonation  point,  on  interaction  with 
the  air-water  interface  directly  below  the  detonation,  excite  the  usual 
plane  waves  in  the  water,  which  are  then  propagated  at  all  angles  of 
incidence  (b)  satisfying  the  condition  cosB>l/n,  where  n  is  the  index  of 
refraction  at  the  interface.  It  is  only  generated  when  the  altitude  of  the 
detonation  is  such  that  at  the  ocean  surface  directly  below  the  detonation 
the  shock  wave  velocity  is  greater  than  the  speed  of  sound  in  the  water 
( n >1 ) .  In  the  design  of  the  classification  system  it  is  assumed  that  the 
contribution  of  the  lateral  wave  energy  to  the  total  signal  received  at  long 
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range  is  negligible.  That  is,  only  the  energy  arriving  via  ray  acoustic 
paths  is  treated.  The  ray  acoustic  paths  of  interest  are  those  associated 
with  the  long  range  propagation  modes.  These  are  illustrated  in  Figure  1. 

The  angle  of  incidence  (e.)  at  the  air-water  interface  is  related  to  the 
angle  of  transmission  (0^.)  at  the  interface  by 

0.  =  arccos  (ncosQ-t)  (4) 

where  n  is  the  index  of  refraction.  The  value  of  n  at  the  interface  is 
given  by 

incident  shock  wave  velocity 

speed  of  sound  in  water  at  the  surface  ' 

Since  the  shock  wave  velocity  varies  (decreases)  with  range  from  the 
detonation  point,  the  value  of  n  at  the  interface  varies  (decreases)  with 
slant  ranqe  to  the  interface,  e.,  .  ,  is  a  function  of  the  near-surface 
sound  speed  maximun  and  is  normally  very  small.  When  the  maximun  occurs  at 
the  surface  (no  surface  duct),  9t(m1n)  is  0°.  et(max)  1S  a  function  of 
sound  speed  at  the  ocean  bottom  and,  for  deep-water  conditions,  is  typically 
on  the  order  of  6  or  7  degrees. 

The  first  computational  step  in  the  determination  of  the  yield-vs- 
altitude  curve  is  the  development  of  a  family  of  relative  received  level 
curves.  Each  of  these  curves  shows,  for  a  constant  yield,  received  signal 
level  as  a  function  of  altitude  relative  to  that  which  would  be  received  if 
the  source  were  located  at  some  reference  position  on  the  depth/altitude 
axis.  The  reference  position  as  defined  here  is  that  depth  within  the  water 
column  at  which  the  propagation  loss  to  the  receiver  location  is  minimized. 
For  a  source  in  deep  water,  this  depth  will  be  that  of  the  deep  sound  channel 
axis.  An  equation  for  determining  relative  received  level  for  the  deepwater 
case  was  developed  in  Reference  1  and  is  presented  below. 


11 


Detonation  Point 


Sound  Speed 


Water 


-  RSR  - 
Region 


Critical  Depth 


Bottom 
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RRL  =  10  log  ^9.255  x  10"13[n2/ ( 1-n2  )  ][P s2 a 2W_0 * 75 3  ]|>s3 /©S(;] 

where, 

RRL  is  the  relative  received  level  in  dB. 
s  is  the  slant  range,  in  feet,  from  the  air  detonation  to  the  region  on 
the  ocean  surface  where  the  ray  acoustics  energy  enters  the  water, 
n  is  the  index  of  refraction  at  the  air-water  interface  at  slant  range  5. 
P$  is  the  peak  overpressure,  in  psi ,  at  slant  range  5. 

W  is  the  detonation  yield  in  kilotons. 

e$  is  the  maximum  transmission  angle  at  the  surface  for  RSR 
propagation,  in  radians.  ,, 

e_  is  the  maximum  angle  of  inclination  at  the  sound  channel  axis 
associated  with  the  long  range  propagation  modes  (RRR  and  RSR),  in 
radians. 

Eq. (6)  is  based  on  the  assumption  that  the  index  of  refraction  does  not 

vary  with  range  across  the  region  within  which  the  energy  propagated  via  RSR 

paths  enters  the  water.  At  high  altitudes,  where  n«T ,  the  assumption  of  a 

/( 

constant  n  is  valid.  However,  for  near-surface  detonations,  as  n  approaches 
1,  the  equation  breaks  down  and  corrections  to  account  for  variations  in  n 
must  be  made. 

In  Reference  1,  numerical  values  of  RRL  based  on  Eq.  (6)  were  provided, 
assuming  a  standard  deep-ocean  environment.*  The  assumed  sound  speed  as  a 
function  of  depth  is  shown  in  Figure  2.  The  sound  speed  is  5037  ft/sec  at  the 
surface,  4900  ft/sec  at  the  sound  channel  axis  (4000  ft)  and  5074  ft/sec  at 

*It  should  be  noted  that  the  derivation  of  Eq.  (6)  involved  the  assumption  of 
a  standard  atmosphere,  the  properties  of  which  are  described  in  Reference  1. 
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the  bottom  (16,400  ft).  The  depth  excess  is  approximately  2640  ft.  The 

values  of  ©s  and  ©sc  for  this  profile  are  6.92°  and  15.05°,  respectively. 

Figure  3  shows  the  results  for  a  1-kiloton  detonation.  The  curve  has  been 

corrected  at  the  lower  altitudes  to  account  for  variations  in  n  with  range. 

1/3 

Results  for  any  other  yield  may  be  obtained  by  a  simple  W  scaling 

procedure.  For  a  yield  W,  the  altitude  of  each  point  on  the  curve  is 
1/3 

multiplied  by  W  .  Then,  each  point  on  the  new  curve  is  shifted  in  level 

by  adding  -0.87  log  W  dB  to  the  relative  received  level. 

The  family  of  RRL  curves  will  consist  of  thirty-six  curves  extending  from 
-2  5 

10  to  10  kilo  tons.  The  yield  associated  with  each  curve  will  differ  by 
a  factor  of  10^5  from  the  yields  for  adjacent  curves.  Thus,  in  ascending 
order  beginning  with  0.01  kilotons,  the  yields,  in  kilotons,  associated  with 
the  first  eleven  curves  are:  0.01,  0.0158  ,  0.0251,  0.0398  ,  0.0631  ,  0.1,  0.158  , 
0.251,  0.398,  0.631  and  1.  Using  the  standard  ocean  environment,  the  RRL 
curves  for  every  fifth  yield  in  the  family  have  been  determined  (scaled  from 
the  1-kiloton  curve)  and  are  presented  in  Figure  4.  Each  curve  is  terminated 
at  an  altitude  of  5000  feet,  which  corresponds  to  the  approximate  maximum 
altitude  of  applicability  of  the  standard  atmosphere  assuned  in  the 
theoretical  developments  upon  which  this  study  is  based. 

From  Eq.  (6)  it  is  seen  that  the  RRL  curve  varies  as  the  cube  of  e$  _  a 
strong  functional  dependence.  Since  ©s  is  a  sensitive  function  of  surface 
sound  speed  and  bottom  depth  and  can  exhibit  marked  spatial  and  temporal 
variability,  it  is  evident  that  large  variations  from  the  results  given  in 
Figure  4  can  occur,  due  to  variations  in  this  one  parameter,  depending  on  the 
time  and  location  of  a  detonation.  As  a  conservative  example,  in  one 
deep-water  region  off  the  U.S.  east  coast,  it  is  estimated  that  over  the 


course  of  a  year  @s  varies  from  4.26°  to  11.3°.  This  corresponds  to 
variations  in  RRL  of  13  dB. 

On  a  worldwide  basis,  surface  sound  speeds  vary  by  no  more  than  about  7% 
throughout  the  year.  Hence,  the  value  of  n  for  a  given  peak  overpressure  will 
never  vary  by  more  than  about  7%.  The  values  of  6  and,  hence,  P  ,  which 
depend  on  surface  sound  speed,  also  are  not  expected  to  change  appreciably,  on 
a  percentage  basis,  from  one  environment  to  another.  Thus,  the  expectation  is 
that  for  any  environment  other  than  that  for  the  standard  ocean,  the 
difference  in  RRL  for  any  given  altitude  will  be  determined  almost  entirely  by 
the  change  in  from  that  for  the  standard  ocean,  with  a  relutively  small, 
but,  perhaps  significant  contribution  due  to  the  change  in  es(;.  And, 
indeed,  sample  calculations  performed  by  the  authors  demonstrate  this  to  be 
the  case.  This  permits  a  simple  scaling  procedure  to  be  employed  whereby  the 
results  obtained  for  the  standard  ocean  can  be  converted  to  RRL  for  any  other 
environment  simply  by  adding  a  correction  value,  involving  9$  and  esc 
only,  to  the  standard  ocean  levels.  This  correction  is  given  by 

A  RRL  =  21.73  +  10  log  (©|/©sc)  (7) 

where 

aRRL  is  the  correction  factor  in  dB 

es  and  esc  are  the  surface  and  sound  channel  angles  for  the  new 

environment,  in  radians 

Since  aRRL  is  independent  of  altitude  the  effect  of  the  scaling  is  simply 
the  shift  of  the  entire  standard  ocean  relative  received  level  curve  by  aRRL. 

The  methodology  to  be  employed  for  the  determination  of  e$  and  esc 
for  a  given  event,  described  below,  is  based  on  the  assumption  that  the 
geometric  spreading  of  the  shot -generated  acoustic  energy  in  the  vicinity  of 
the  receiving  sensor  is  essentially  cylindrical;  that  is,  that  the  receiver  is 

1 
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at  long  range  from  the  detonation.  This  assumption  is  not  believed  to  detract 
from  the  general  usefulness  of  the  systan  in  view  of  the  strong  expectation 
that,  for  a  given  event,  the  detonation/receiver  separation  range  will  be  very 
large  -  perhaps  thousands  of  miles. 

Long  range  propagation  loss  within  the  deep  sound  channel  may  be 
approximated  by  assuming  spherical  spreading  out  to  some  transition  range, 
followed  by  cylindrical  spreading  thereafter.  The  loss,  as  a  function  of 
range,  is  given  by 

PL  =  10  log  Rq  +  10  Log  R  +  aR  (8) 

where 

PL  is  the  propagation  loss  in  dB 
Rq  is  the  transition  range  in  yards 
R  is  the  source/receiver  separation  range  in  yards 
a  is  the  absorption  coefficient  in  dB/yard 
The  transition  to  cylindrical  spreading  results  from  the  interaction  of  the 
acoustic  energy  propagating  outward  from  the  source  with  the  ocean  surface  and 
bottom  boundaries  which  serve  to  trap  the  energy  within  the  water  column.  In 
reality,  of  course,  the  transition  from  spherical  to  cylindrical  spreading  is 
gradual  and  does  not  occur  abruptly  at  a  discrete  range.  RQ  is  simply  a 
range  within  the  transition  region  which  yields  reasonable  predictions  at  long 
range.  Consider  the  illustration  shown  as  Figure  5.  In  Figure  5  the  source 
is  located  at  mid-column  for  illustrative  purposes  only.  The  angle  © 
represents  the  maximum  inclination  angle  associated  with  the  long  range 
propagation  modes  as  would  be  determined  from  the  sound  speed  profi  le.  The 
dashed  curves  represent  the  refracted  paths  associated  with  ©,  while  the 
straight  line  puths  ignore  refraction.  For  this  study,  RQ  ,s  taken  as  that 
range  at  which  the  vertical  separation  of  the  straight  line  paths  associated 
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Ro  H  Ocean  Surface 


Figure  5.  Diagram  Illustrating  the  Relationship  Between  the 
Spherical  Spreading-to-Cylindrical  Spreading 
Transition  Range  (Rq) ,  the  Water  Depth  (d)  and  the 
Maximum  Inclination°Angle  at  the  Source  (0)  Associated 
With  the  Long-Range  Propagation  Modes.  (See  text.) 


with  e  equals  the  water  depth  (d).  This  reasonably  estimates  the  range  at 
which  the  channeling  effect  of  the  ocean  boundaries  begins  to  take  effect. 

Rq  is  then  given  by 

Rq  =  d/2  tan©-  (9) 

Substituting  the  expression  for  Rfl  in  Eq.  (8),  we  obtain 

PL  =  10  log  (d/2tane)  +  10  log  R  ♦  aR  (10) 

For  a  known  propagation  loss  obtained  at  long  range,  an  estimate  of  the  angle 
©  at  the  source  may  be  obtained  by  solving  Eq.  (10)  for  9.  The  resulting 
equation  is  given  by 


If  the  known  propagation  loss  is  obtained  for  a  surface  source,  a  3-dB 
correction  must  be  applied  to  Eq.  (10)  to  account  for  the  absence  of  upgoirig 
energy  at  the  source.  The  loss  is  then  given  by 

PL  =  10  log  (d/2tane)  +  10  log  R  +  aR  +  3  (12) 

and  9  is  given  by 


The  values  of  9$(.  and  9S  for  use  in  the  determination  of  the  RRL  curves 
for  the  atmospheric  case  will  be  determined  using  Eqs.  (11)  and  (13), 
respectively.  The  PL  inputs  to  the  equations  are  obtained  from  the  set  of 
propagation  loss  computations  which  were  obtained  for  assumed  source  positions 
throughout  the  water  column  (Section  2.1.1).  The  PL  input  for  the 
determination  of  ©sc  (Eq.  (ll))  is  the  loss  obtained  for  an  on-axis  source  - 
if  a  deep  channel  axis  is  present.  If  a  deep  axis  is  not  present,  the  PL 
input  to  Eq.  (11)  is  the  minimum  predicted  loss  within  the  water  column.  The 
value  of  PL  to  be  used  in  the  determination  of  ©s  is  that  obtained  for  a 


surface  source.  The  water  depth  (d)  is  assumed  to  be  that  at  the  receiver 
site.  The  values  obtained  for  esc  and  &s  are  equivalent  to  the  angles 
which  would  be  obtained  for  the  predicted  PL  values  if  the  water  depth  along 
the  entire  propagation  track  (detonation  site  to  receiver)  was  equal  to  that 
at  the  receiver  site  and  if  the  acoustic  bottom  loss  along  the  track  was 
infinite.  The  characteristics  of  the  actual  environment  (sound  speed 
structure,  bottom  loss,  bathymetry,  etc.)  are  reflected  in  the  computed  values 
of  o$  and  ©sc  through  the  PL  input. 

The  determination  of  the  RRI.  curves,  then,  comprises  three  steps.  Step  1 

is  to  scale  the  1-kiloton  standard  ocean  RRL  curve  to  the  thirty-five  other 

1/3 

yields  in  the  RRL  curve  family  (scaling:  altitude  multiplied  by  W 
followed  by  a  shift  of  -O.871og  W  dB) .  Step  2  is  to  determine  e  and  &sc 
for  the  subject  environment  using  Eqs.  (11)  and  (13).  Step  3  is  to  scale  the 
standard  ocean  RRL  curves  obtained  in  Step  1  to  the  subject  environment,  using 
Eq.  (7). 

The  next  computational  step,  once  the  RRL  curves  in  the  subject 
environment  have  been  determined,  is  to  normalize  the  RRL  curves  such  that, 
for  any  given  normalized  curve  (NRRL),  the  NRRL  values  represent  levels 
relative  to  the  level  which  would  be  received  for  a  1-kiloton  detonation  on 
axis  (or,  the  depth  at  which  the  detonation  site-to-recei ver  propagation  loss 
is  minimized).  For  each  RRL  curve,  the  normalization  is  performed  by  adding 

7.5  log  W  dB  to  the  RRL  values,  where  W  is  the  yield,  in  kilotons,  for  which  a 

given  curve  was  computed.  This  results  in  a  shift  of  7.5  log  W  dB  to  the 

entire  RRL  curve.  For  example,  for  the  100-kiloton  RRL  curve,  the 

normalization  factor  is  7.5  log  100,  or,  +15  dB.  At  any  given  altitude,  then, 
NRRL  =  RRL  ♦  15  dB.  The  correction  for  each  curve  merely  reflects  the 
difference  in  source  level  between  that  of  an  underwater  shot  of  yield  W  and 


that  of  an  underwater  1-kiloton  detonation.  Corrections  range  from  -15  dB  for 
the  0.01-kiloton  curve  to  +37  .5  dB  for  the  10^-kiloton  curve. 

The  determination  of  the  yield-vs-al  t  itude  curve  for  the  subject 
environment  is  based  directly  on  the  NRRL  curve  family.  Let  us  assume,  for 
discussion  purposes,  that,  for  a  particular  case  the  subject  environment  is 
identical  to  that  of  the  standard  ocean  environment.  The  NRRL  curves  for  this 
case  are  as  shown  in  Figure  6.  For  presentation  purposes,  only  every  fifth 
curve  in  the  total  family  of  curves  is  presented  in  the  figure.  Let  us  also 
assure  that  on  the  underwater  yield-vs-depth  curve  for  a  particular  event 
occurring  in  this  environment,  the  yield  at  the  reference  depth  has  been 

Q 

determined  to  be  0.02  lb  (10  kiloton).  The  received  level  for  this  yield 
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is,  then,  -60  dB  (7.5  log  10  )  relative  to  the  level  that  would  be  received 

if  the  yield  at  the  reference  depth  were  1  kiloton.  Then,  for  any  given  NRRL 
curve  the  yield/altitude  pairing  at  the  point  on  the  curve  for  which  NRRL 
equals  -60  dB  defines  a  point  on  the  yield-vs-alt  itude  curve,  and  the  set  of 
points  defines  the  yie ld-vs-al  titude  curve.  In  Figure  6,  the  yie ld-vs-al  titude 
curve  for  this  example  is  determined  by  the  points  indicated  by  the  heavy  dots 
along  the  vertical  line  at  -60  dB.  The  yie  ld-vs-al  titude  plot  for  this  data 
is  shown  in  Figure  7.  The  curve  definition  is,  of  course,  not  nearly  as  good 
as  it  vould  be  if  the  full  set  of  NRRL  curves  was  used. 

It  is  important  to  note  that  the  1-kiloton  normalization  may  be  performed 
before,  rather  than  after,  the  environmental  correction  to  the  standard  ocean 
(Eq.  (7))  has  been  applied.  That  is,  the  normalization  may  be  performed  on 
the  standard  ocean  RRL  curve  family,  and  the  environmental  correction  applied 
to  the  resulting  standard  ocean  NRRL  curves  to  obtain  NRRL  in  the  subject 
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environment.  When  the  classification  system  is  implemented,  the  NRRL  curves 
for  the  standard  ocean  will  be  permanently  stored  within  the  computer.  Then, 
for  a  given  application  of  the  system  the  only  computational  tasks  required 
for  the  determination  of  the  yield-vs-alt  itude  curve  in  the  atmosphere  are: 

(1)  determine  es  and  &sc  for  the  subject  environment;  (2)  apply  the 
environmental  correction,  Eq.  (7),  to  the  standard  ocean  NRRL  curves  to  obtain 
NRRL  curves  for  the  subject  environment;  (3)  determine  the  received  level  for 
the  reference-depth  yield  on  the  underwater  yield-vs-depth  curve  relative  to 
that  which  would  be  received  by  a  1-kiloton  detonation  at  the  reference  depth; 
(4)  specify  the  yield-vs-depth  curve  by  the  determination  of  yield/ altitude 
points  within  the  NRRL  curve  family  for  which  the  value  of  NRRL  is  that 
determined  in  the  previous  step. 

2.1.3  Yield  Vs  Depth/Altitude  -  Near-Surface  Case 
2. 1.3.1  Above-Surface  Region 

In  this  section  we  specify  the  methodology  for  the  determination  of  yield 
vs  altitude  over  the  region  of  the  altitude  axis  for  which  the  height 
associated  with  any  given  yiald  on  the  yie ld-vs-al  titude  curve  is  less  than 
the  range  for  hydrodynamic  separation.  Within  the  near-surface  region  the 
water  participates  in  the  formation  of  the  fireball  and  a  shock  wave  is 
generated  in  the  water  medium  as  in  the  case  of  an  underwater  detonaton.  This 
is  a  very  efficient  mode  of  excitation  since  it  does  not  involve  the  losses 
associated  with  the  impedance  mismatch  between  air  and  water.  A  complete 
discussion  of  the  theory  underlying  the  methodology  to  be  employed  here  will 
not  be  presented.  Rather,  the  reader  is  referred  to  Reference  3  which 
presents  a  complete  treatment. 
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The  first  step  in  the  computational  process  is  the  development  of  a 

family  of  RRL  curves.  As  for  the  atmospheric  case,  the  RRL  curve  family  will 

-2  5 

consist  of  thirty-six  curves  extending  from  10  to  10  kilotons,  with 
each  curve  differing  by  a  factor  of  10^  in  yield  from  the  yields  for 
adjacent  curves. 

For  a  shot  occurring  in  the  air  at  an  altitude  less  than  the  range  for 
hydrodynamic  separation,  there  is  penetration  of  the  fireball  into  the  water. 
The  depth  of  penetration  is  inversely  related  to  shot  altitude.  The  fraction 
of  the  total  energy  available  for  formation  of  the  shock  wave  in  the  water  is 
just  that  fraction  of  the  fireball  energy  which  penetrates  the  air-water 
interface.  This  directly  induced  shock  wave  energy,  normalized  to  the  induced 
energy  of  a  detonation  occurring  at  the  air-water  interface,  is  shown  as  a 
function  of  scaled  detonation  altitude  in  column  2  of  Table  1.  The  scaled 
altitude  values  are  equal  to  altitude  (H)  divided  by  hydrodynamic  separation 
range  (R).  The  values  range  from  0  dB  at  the  normalization  point  to  -oodB 

a 

at  the  altitude  corresponding  to  Ra,  and  are  assumed  to  be  environmentally 
independent.  The  RRL  for  a  surface  detonation  is  assumed  here  to  be  -11.1  dB, 
independent  of  yield  and  environmental  conditions.  The  determination  of  this 
value  is  based  on  three  considerations.  First,  using  the  standard  ocean 
environment,  it  was  determined  that,  based  on  the  solid  angles  at  the  source 
associated  with  long  range  propagation,  the  received  level  for  an  explosive 
source  just  below  the  ocean  surface  is  -3.4  dB  relative  to  the  level  that 
would  be  received  from  an  on-axis  source  of  equal  yield.  For  a  source  right 
at  the  surface  the  received  level  will  be  reduced  an  additional  3  dB  due  to 
the  absence  of  initially  upgoing  energy  at  the  source,  yielding  a  relative 
level  of  -6.4  dB.  Based  on  data  obtained  for  small  HE  charges  at  the 
air-water  interface  it  is  estimated  that  directly  below  a  surface  charge  the 
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Table  1.  Directly  Induced  Shock  Wave  Energy  as  a 
Function  of  Scaled  Detonation  Altitude 


I 

I 

I 

I 

I 

I 

I 

I 

I 


(1) 

Scaled  Detonation 
Altitude  (H/Ra) 

0 

0.1 


0.2 


0.3 

0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1.0 


(2)  (3) 

Directly  Induced  Shock  Wave 
Energy  Relative  to  That  of  a 
Surface  Detonation  (dB)  RRL  (dB) 

0.0  -11.1 

-  0.7  -11.8 

-  1.5  -12.6 

-  2.5  -13.6 

-  3.7  -14.8 

-  4.9  -16.0 

-  6.4  -17.5 

-  9.2  -20.3 

-12.5  -23.6 

-18.4  -29.5 

-  oo  -  oo 


i: 


shock  wave  energy  flux  will  be  34%  of  the  t‘norgy  for  an  underwater  detonation. 

Thus,  an  additional  correction  of  -4.7  dB  is  required,  yielding  the  total  RRL 

of  -11.1  dB  tor  a  surface  detonation.  It  is  seen  from  the  above  discussion 

that  RRL  for  a  surface  source  is  not  truly  environmentally  independent. 

However,  it  is  not  expected  that  this  assumption  will  lead  to  significant 

errors  in  foreseeable  real-world  situations. 

RRL  as  a  function  of  scaled  altitude  is  shown  in  column  3  of  Table  1. 

These  values  were  obtained  simply  by  adding  -11.1  dB  to  the  relative  directly 

induced  energy  levels  given  in  column  2.  RRL  as  a  function  of  absolute 

altitude  for  a  1-kiloton  detonation  is  shown  in  Figure  8.  These  results  may 

1/3 

be  scaled  to  any  other  yield  by  simple  W  '  scaling.  That  is,  for  any  yield 

W,  expressed  in  kilotons,  the  altitude  of  each  point  on  the  1-kiloton  curve  is 
1/3 

multiplied  by  W  .  The  RRL  curves  for  every  fifth  yield  in  the  curve 
family  are  presented  in  Figure  9.  It  should  be  noted  that  each  curve  is 
terminated  at  that  altitude  for  which  the  value  of  RRL  is  equal  to  RRL  at  the 
lowest  altitude  (hydrodynamic  separation  range)  on  the  corresponding 
atmospheric  RRL  curve  (see  Figure  4).  The  implemented  classification  system 
will  simply  ignore  the  small  gap  introduced  by  this  truncation  as  its 
inclusion  can  produce  spurious,  physically  unmeaningful  results  in  the 
determination  of  yield  vs  altitude. 

Once  the  RRL  curve  family  has  been  developed,  the  process  of  determining 
yield  vs  altitude  in  the  near-surface  region  proceeds  in  precisely  the  same 
fashion  as  for  the  atmospheric  case  except  that  there  is  no  environmental 
correction  to  consider.  Specifically,  the  remaining  computational  steps  are: 
(1)  develop  a  family  of  NRRL  curves;  (2)  determine  the  value  of  NRRL  for  the 
reference  depth  yield,  and  (3)  specify  the  yiel d-vs-al ti tude  curve  in  the 
near-surface  region  by  the  determination  of  yield/altitude  points  within  the 
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NRRL  curve  family  for  which  the  value  of  NRRL  is  that  determined  for  the 
reference  depth  yield.  NRRL  has  precisely  the  same  definition  here  as  given 
in  Section  2.1.2  and  is  determined  in  precisely  the  same  fashion.  NRRL  curves 
for  the  near-surface  region  are  presented  in  Figure  10.  Again,  for 
presentation  purposes,  only  every  fifth  curve  in  the  family  of  curves  is 
presented . 

Sample  yield-vs-al titude  curves  combining  the  near-surface  and 
atmospheric  regions  are  presented  in  Figure  11.  For  the  atmospheric  region, 
determinations  were  based  on  the  standard  environment.  Curves  are  shown  for 
seven  reference  depth  yields  ranging  from  0.02  lb  to  1  kiloton.  The  curves 
are  based  on  the  set  of  NRRL  curves  consisting  of  every  fifth  curve  of  the 
near-surface  and  atmospheric  curve  families.  When  the  classification  system 
is  implemented,  tne  near-surface  NRRL  curves  will  be  permanently  stored  within 
the  computer.  Then,  for  a  given  application  of  the  system  the  only 
computational  tasks  required  for  the  determination  of  yiel d-vs-al titude  in  the 
near-surface  region  will  be  steps  2  and  3  above. 

The  development  of  the  methodology  for  determination  of  yield  vs  depth/ 
altitude  has  led  to  an  important  general  conclusion  relative  to  the 
classification  of  a  detonation.  If  we  consider  the  practical  minimum  yield  of 
a  test  nuclear  charge  to  be  about  10  tons  (NRRL  =  -15  dB),  then,  for  a  given 
detonation-like  signal  which  indicates  a  possible  underwater  nuclear 
detonation  (>10  tons  on  the  y iel d-vs-depth  curve),  one  can  conclude  with  a 
high  degree  of  confidence  that  the  signal  was  not  generated  by  an  atmospheric 
blast  (altitude  >15W^/^  ft).  The  rationale  for  this  conclusion  is  that  the 
generation  of  such  a  large  underwater  signal  would  require  an  impractically 
high  atmospheric  yield.  However,  the  possibility  that  the  detonation  occurred 
in  the  near-surface  region  could  not  be  ruled  out. 
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2. 1.3. 2  Belo w-Surf ace  Region 

As  previously  discussed,  the  underwater  near- surface  region  is  defined  as 

1/3 

that  portion  of  the  water  column  shallower  than  21  W  feet,  where  W  is  the 
detonation  yield  in  kilotons.  Within  this  region  all  of  the  water  above  the 
detonation  is  vaporized.  To  date,  no  phenomenologically  based  methodology  for 
the  determination  of  RRL  at  near- surface  depths  has  been  devised.  The 
classification  system  will  determine  yield  vs  depth  within  the  near- surface 
region  by  simple  linear  interpolation.  That  is,  yield/depth  pairings  will  be 
linearly  interpolated  based  on  the  yield/depth  combination  at  the  terminus  of 
the  underwater  yield-vs-depth  curve  and  the  yield  obtained  for  a  surface 
source.  The  surface  yield  is  given  by 

WSIR  =  30,2  WRD 

where 

WSUR  the  surface  yield 
WRp  is  the  reference  depth  yield 

2.2  Signal  Processing 

The  signal  from  a  nuclear  underwa  ter/ ai  r  burst  is  assumed  to  be  received 
on  a  single  hydrophone  located  within  the  deep  sound  channel.  It  is  further 
assumed  that  the  underwater  sound  path  between  the  burst  and  the  receiving 
hydrophone  is  long  (~2000  nm)  and  will  support  RRR  and/or  RSR  propagation. 

The  time  constant  of  the  positive  overpressure  pulse  at  the  receiving 
hydrophone  is  short  (~one  second,  maximum),  but,  because  of  the  arrival  time 
dispersion  of  the  signal  propagating  along  the  many  multi  paths ,  the  received 
signal  consists  of  a  series  of  overlapping  pulses  that  at  long  ranges  may 
persist  for  many  seconds.  This  time  dispersion  varies  with  range;  at  a  range 
of  2000  nm  in  the  deep  ocean,  it  is  approximately  15  seconds.  The  operational 
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frequency  passband  of  the  receiver  is  10  to  30  Hz,  a  bandwidth  of  about  1.6 
octaves.  The  signal  within  this  passband  will  have  a  spectrin  slope  of 
-6  dB  octave  for  underwater  bursts  and  for  high  yield  air  bursts,  a  spectrum 
slope  of  0  dB/octave  for  moderate  yield,  moderate  altitude  air  bursts  and  a 
spectrum  slope  of  +6  dB/octave  for  low  yield  air  bursts.  Most  of  the  yield/ 
altitude  cart inations  for  which  bursts  can  be  detected  fall  into  the  S  dB/octave 
spectrum  slope  category  which  will  be  the  basis  of  the  following  signal 
processing  analysis.  The  absorption  loss  in  the  10  to  30  Hz  passband  is 
negligibly  small  and  has  no  effect  on  the  signal  spectrum  slope. 

The  initial  signal  detection  can  best  be  obtained  by  a  broadband  square 
law  (auto-correlator)  detector.  The  output  signal-to-noi  se  likelihood  ratio 
is  given  by 


where 

SiN  is  the  signal  spectrum  level  at  the  detector  input 
Njn  is  the  noise  spectrum  level  at  the  detector  input 
F  is  the  passband  shaping  filter  function 

is  the  lover  frequency  of  the  passband 
^2  is  the  upper  frequency  of  the  passband 
T  is  the  averaging  time  of  the  detected  output  signal 

o 

If  the  shaping  filter  is  proportional  to  SjN/NjN  .  it  is  the  optimum  filter 
and  Eq.  (15)  reduces  to 
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The  deep  water  noise  level  in  the  10  to  30  Hz  passband  for  an  area  of 
li ght-to-moderate  shipping  noise  is  approximately  given  by 

N,n  =  86  -  10  log  f,  dB  re  luPa/Hz  (17) 

The  signal  level  at  the  hydrophone  receiver  for  an  underwater  nuclear 
detonation  at  the  depth  for  minimum  propagation  loss  is 

SIN  =  SL  -  PLmin,  dB  re  luPa/Hz  (18) 

where 

SL  is  the  spectrun  source  level  of  the  detonation  in  dB  re  lyPa/Hz  at 
1  yard 

PLmin  is  the  propagation  loss  from  source  to  hydrophone  receiver  in  dB 
The  spectrum  source  level  (SL)  in  the  10  to  30  Hz  passband  for  a  1-kiloton 
underwater  detonation  is  294  -  20  log  f,  dB.  Since  the  normalized  relative 
received  level  (NRRL)  gives  the  source  level  of  any  nuclear  detonation 
relative  to  an  underwater  detonation  of  1-kiloton  yield,  the  spectrum  source 
level  can  be  expressed  as 

SL  =  294  -  20  log  f  +  NRRL,  dB  re  lyPa/Hz  at  1  yard  (19) 

Thus,  the  signal  level  at  the  receiver  input  is 


SIN  =  294  -  20  log  f  +  NRRL  -  PLMIN,  dB  re  lyPa/Hz;  (20) 

the  si gnal-to-noise  ratio  at  the  receiver  input  is 

($)  =  208  -  10  log  f  +  NRRL  -  PLMIN,  dB,  (21) 

and  the  detected  output  signal-to-noi se  ratio  for  an  averaging  time  T  of  15 

seconds  and  a  passband  of  10  to  30  Hz  is,  from  Eqs.  (16)  and  (21), 


416  -  2  PlMJN  ♦  2  NRRl  +  10  log  15 
416  -  2  PLmin  +  2  NRRL,  dB 


/ 


10 


(22) 


The  minimum  level  of  NRRL  that  can  be  reliably  detected  is  determined  by  a 
desired  probability  of  detection,  together  with  an  allowed  probability  of 
false  alarm  per  averaging  time  interval  (T).  A  measure  of  these  joint 
probabilities  is  a  parameter  called  the  detection  index  (d)  which  can  be 
obtained  from  a  set  of  Receiver  Operating  Characteristic  (ROC)  curves.  The 
conventional  ROC  curves  assume  a  steady  signal  over  the  averaging  time, 
stationary  Gaussian  noise,  a  large  time-bandwidth  processing  gain  and  no 
interfering  signal-like  noises.  Whenever  these  conditions  are  not  all  met,  a 
higher  detection  index  is  required.  For  this  application  (the  detection  of 
nuclear  detonations)  the  signal  is  usually  not  steady  over  the  averaging  time 
because  of  the  many  mul ti paths  with  varying  transmission  losses;  the  noise  in 
the  10  to  30  Hz  passband  is  both  non-s tat ionary  and  too  spiky  to  be  Gaussian; 
the  time-bandvidth  gain  is  too  small;  and,  shipping  noises  will  often  cause 
interference  problems.  Because  of  all  these  factors,  a  rather  large 

detection  index  is  recommended.  Using  Eq.  (22),  let  =  d  =  23  dB  and 

v  "'OUT 

solve  for  the  minimun  detectable  value  of  NRRL.  Thus, 

NRRLMIN  =  7  <23  "  416  +  2  PLMIN)  =  PLMIN  ~  196'5’  dB  (23) 

The  absolute  value  of  signal  output,  using  Eq.  (20)  for  -js 

f2 

S0UT  =  20  lo9  T /  BIN  df>  dB  re  lwPa2 
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I#:  r\ 

=  588  -  2  PLmin  +  2  NRRL  +  20  log  15  f  dB  re  lyPa^ 

10  f 

=588-2  PLmin  +  2  NRRL,  dB  re  l^Pa2  (24) 

At  the  minimun  detectable  level  of  NRRL, 

(S0Ut)min  =  588  -  2  •  196.5  =  195  dB  re  luPa2  (25) 

Whenever  there  is  no  signal  the  detected  and  averaged  output  is  that  due  to 
only  the  ambient  noise.  Its  level  is 

f2  30 

Nout  =  20  log  T j  N1N  df  =  172  +  20  log  15  /  ~  =  189  dB  re  IpPa2  (26) 


Note  that  Nqut  is  6  dB  less  than  (50yj)MIN*  For  a  processor  using  the 
output  of  a  single  am  idirecti  onal  hydrophone  located  in  an  area  where  there 
may  be  many  interfering  noise  sources,  such  as  from  moving  ships,  it  is 

probably  better  to  determine  NRRL  frcm  Sq^^.  and  NQUT  than  from  .  Thus, 

from  Eqs.  (24)  and  (25),  we  obtain 

NRRL  =  (SQUT  -  NqUT)  +  PLmin  -  199.5,  dB  (27) 

In  the  standard  ocean  the  long  range  minimum  propagation  loss  is 
approximately  given  by 

PLmin  =  10  log  5.8  nm  +  10  log  R  nm  +  20  log  2000  yd/nm 

=  73.7  +  10  log  R  nm,  dB  (28) 

At  R  =  2000  nm,  PLMjn  =  106.7  dB.  Therefore,  at  a  range  of  2000  nm  in  the 
standard  ocean,  the  minimim  detectable  value  of  NRRL  using  Eqs.  (25),  (26)  and 
(27)  is  equal  to  -89.8  dB. 

As  an  illustrative  example,  assume  that  a  signal  has  been  received  in  the 
standard  ocean  environment.  Let  the  output  level  of  the  signal,  Sgyy 
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87  d B  above  the  output  level  of  the  ambient  noise  before  the  signal  arrives, 
and  let  the  minimun  propagation  loss  at  a  range  of  1700  nm  be  106  dB .  Then, 
solving  £q.  (27),  IWRL  =  -50  dB.  This  value  of  NRRL  can  now  be  used  in 
Figured  to  determine  possible  values  of  yiel  d/ al  t  itude  for  an  air  nuclear 
detonation.  For  example,  a  1-kilo  ton  yield  would  be  at  an  altitude  of  about 
100  ft,  a  10-kilton  yield  at  about  450  ft,  and  a  100-kiloton  yield  at  about 
1500  ft.  This  could  not  have  been  an  underwater  nuclear  detonation  because 
NRRL  is  too  smal 1 . 

After  the  detection  of  a  nuclear  detonation  has  been  made,  it  might  be 
desirable  to  perform  a  spectral  analysis  of  the  received  signal  in  order  to 
determine  whether  the  higher  critical  frequency  is  below,  within  or  above  the 
10  to  30  Hz  passband.  A  1-Hz  narrowband  resolution  and  an  averaging  time  of 
15  seconds  are  probably  suitable  for  obtaining  the  signal  spectrum.  If  this 
spectrun  has  a  slope  of  -6  to  -3  dB/octave,  the  critical  frequency  is  below 
the  passband  and  no  correction  to  NRRL  due  to  spectrum  shape  is  required.  If 
this  spectrun  has  a  slope  of  -3  to  +6  dB/octave,  NRRL  is  too  low  and  should  be 
corrected  to  a  higher  value.  The  amount  of  this  correction  is  not  large  if 
the  spectrun  slope  is  less  than  0  dB/octave,  but  if  it  is  greater  than  0  dB/ 
octave  the  correction  is  likely  to  be  significant.  However,  only  low  yield, 
low  altitude  nuclear  detonations  in  the  air  will  have  a  positive  spectrun 
slope.  This  is  the  yield/  al  titude  region  which  is  subject  to  large  NRRL 
errors  due  to  other  effects.  It  is  also  unlikely  that  such  low  yield  air 
detonations  will  be  tested  and  detected.  Therefore,  it  is  probably  not 
necessary  to  provide  corrections  to  NRRL  due  to  spectrun  shape. 


2 . 3  Err  or  Anal  ysis 

The  passband  filtered,  detected  and  averaged  output  of  the 
autocorrelation  processor  for  the  receiver  hydrophone  provides  the  basic 
measurements  that  are  used  to  obtain  yield-vs-depth/ al  titude  estimates.  These 
measurements  are  the  processor  output  level,  in  dB,  due  to  the  received 
acoustic  energy  from  a  nuclear  detonation  (Sq^)  and  the  output  level  due  to 
ambient  noise  just  before  the  arrival  of  the  signal  (Nqut).  The  dB 
difference  in  these  tw>  measurements  (Sgy j  -  Nq^j)  can  be  obtained  quite 
accurately,  provided  the  hydrophone  and  processor  are  linear  over  the  dynamic 
range  of  the  signal.  This  dynamic  range  should  be  of  the  order  of  60  dB  for 
air  bursts  or  underwater  detonations,  and  about  120  dB  for  both.  Such  a  large 
dynamic  range  requi  res  an  accurate  automatic  gain  control  (AGC)  for  the 
hydrophone  receiver  output,  from  which  the  full  dynamic  range  of  the  signal 
can  be  linearly  reconstructed.  The  standard  deviation  error  in  the 
measurement  of  (SQyT  -  NqUT)  is  approximately  3.5  dB  if  the  receiver 
system  is  linear  but  not  calibrated.  About  2.5  dB  of  this  error  may  be  due  to 
an  inaccurate  estimate  of  the  ambient  noise  level.  If  the  receiver  system  is 
calibrated,  this  source  of  error  can  be  eliminated,  which  reduces  the 

estimated  SQUT  measurement  error  to  about  2.5  dB.  With  a  calibrated  system, 
Eq.  (24),  which  does  not  require  a  measurement  of  NqUT  to  determine  NRRL, 
can  be  used.  The  other  source  of  error  in  the  estimation  of  NRRL  is  that  due 
to  error  in  the  propagation  loss  (PL)  calculation,  which  is  estimated  to  be 
about  2  dB  for  ocean  environments  that  approximate  that  of  the  assumed 
standard  ocean.  For  other  types  of  environments,  the  propagation  loss  error 
will  probably  be  greater.  Thus,  the  standard  deviation  error  in  NRRL  is 
estimated  to  be  about  3.2  dB  using  Eq.  (24),  and  4  dB  using  Eq.  (27), 
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depending  upon  whether  or  not  the  receiver  system  is  calibrated,  for  ocean 
environments  that  provide  good  long  range  propagation. 

The  calculation  of  NRRL  does  not  provide  a  solution  of  the  yield  and 
depth/al titude  of  a  nuclear  detonation.  It  simply  reduces  the  total  set  of 
yield-depth/ al titude  pairs  to  those  which  produce  the  same  amount  of  acoustic 
energy  as  that  actu  11 y  received  within  the  passband  at  long  range.  For  an 
underwater  nuclear  detonation,  the  yield  can  be  determined  with  fair  accuracy 
even  if  the  depth  of  the  detonation  is  unknown..  This  is  due  to  the  fact  that 
NRRL  for  a  given  yield  varies  a  maximim  of  about  6.5  dB  with  depth.  However, 
for  an  air  burst,  NRRL  can  vary  by  as  much  as  75  dB  as  a  function  of  altitude 
(from  0  to  5000  feet)  for  a  given  yield.  It  was  convenient  to  divide 
yield/ al titude  as  a  function  of  NRRL  into  two  regions  -  that  for  which  the 
detonation  altitudes  exceed  the  hydrodynamic  separation  range  (Figure  6),  and 
that  for  which  the  fireball  interacts  with  the  water  before  hydrodynamic 
separation  (Figure  10).  It  can  be  seen  from  these  figures  that  the  possible 
yield-altitude  pairs  for  a  given  value  of  NRRL  -  a  vertical  line  -  are  simply 
those  which  intersect  this  line.  Then,  if  by  some  other  means,  the  value  of 
altitude  (or  yield)  can  be  obtained,  the  value  of  yield  (or  altitude)  can  also 
be  obtained  at  the  point  on  the  vertical  line  where  it  intersects  the  altitude 
(or  yield)  line.  An  exception  to  this  statement  can  be  seen  in  Figure  6, 
where,  for  some  values  of  NRRL  and  yield,  the  altitude  is  multivalued.  For 
example,  at  NRRL  =  -50  dB  and  yield  (W)  =  1  kiloton,  the  altitude  may  be 
either  15  feet  or  90  feet. 

Now,  if  the  error  statistics  associated  with  NRRL  are  known,  confidence 
limits  on  yield-vs-altitude  determination  can  be  specified.  Values  within 
these  limits  will  fall  within  an  area  about  the  vertical  line  at  the  mean 
value  of  NRRL.  If  the  NRRL  error  is  a  Gaussian  distributed  random  variable 


with  a  standard  deviation  of  4  dB,  the  probability  that  the  true  yield- 
altitude  pairings  fall  within  the  vertical  lines  at  4  dB  on  either  side  of  the 
mean  value  of  NRRL  is  0.68,  and  the  probability  that  they  fall  outside  these 
lines  is  0.32.  Thus,  if  the  value  of  altitude  (or  yield)  can  be  determined  by 
some  other  means,  the  value  of  yield  (or  altitude)  will  have  a  random  error 
that  is  a  non-linear  function  of  both  the  NRRL  error  and  the  altitude  (or 
yield)  error.  An  inspection  of  Figures  6  and  10  will  show  the  variability  of 
yield  error  as  a  function  of  NRRL  and  altitude  -  called  Zq  in  this 
discussion  -  for  given  values  of  error  in  each.  For  example,  let  NRRL  in 
Figure  6  be  -50  dB  with  a  standard  deviation  error  of  4  dB,  and  let  Zq  be 
1500  feet  with  zero  error.  Then,  there  is  a  50-percent  probability  that  the 
yield  will  be  between  69  and  145  kilotons,  and  a  90-percent  probability  that 
it  will  be  between  40  and  248  kilotons.  (This  assumes  that  the  NRRL  error  is 
a  Gaussian  random  variable  with  zero  mean.)  Next,  let  NRRL  be  -50  dB  with 
zero  error,  and  let  Zq  be  1500  feet  with  a  standard  deviation  error  of  300 
feet.  If  this  altitude  error  is  a  Gaussian  distributed  random  variable  with 
zero  mean,  there  is  a  50-percent  probability  that  the  yield  will  be  between  78 
and  135  kilotons,  and  a  90-percent  probability  that  the  yield  will  be  between 
55  and  208  kilotons. 

Looking  at  another  region  of  yield-altitude  as  a  function  of  NRRL,  let 
NRRL  in  Figure  10  be  -30  dB  with  a  standard  deviation  error  of  4  dB,  and  let 

z0  be  28  feet  with  zero  error.  Then,  there  is  a  50-percent  probability  that 
the  yield  will  be  between  9.7  and  10.4  kilotons,  and  a  90-percent  probability 
that  the  yield  will  be  between  9.2  and  10.9  kilotons.  Now,  let  NRRL  be  -30  dB 

with  zero  error,  and  let  Zq  be  28  feet  with  a  standard  deviation  error  of  6 
feet.  For  this  error  in  altitude  there  is  a  50-percent  probability  that  the 
yield  will  be  between  7.7  and  13.8  kilotons,  and  a  90-percent  probability  that 


the  yield  will  be  between  5.3  and  22.1  kilotons.  These  examples  show  the 
variability  of  the  error  in  yield  as  a  function  of  error  in  NRRL  and  as  a 
function  of  error  in  altitude  for  different  regions  of  al  titude-NRRL  space. 

Another  possible  source  of  error  in  yield-altitude  as  a  function  of  NRRL 
is  the  accuracy  in  the  derivation  and  generation  of  the  curves  in  Figures  6 
and  10.  All  known  variables  were  accounted  for  to  first  order,  but  some  have 
been  extrapolated  to  the  extent  that  second  order  effects  may  introduce 
several  dB  of  error  in  the  acoustic  energy  levels  in  some  regions  of  the 
yield-altitude  space.  Nevertheless,  it  is  believed  that  these  curves  are  as 
accurate  as  can  be  reasonably  expected  from  theory  and  the  available 


3.  SUMMARY 


A  computational  scheme  is  developed  for  the  determination  of  yield  vs 

depth/ al ti tude  based  on  a  received  underwater  detonation-like  acoustic  signal 

in  the  deep  ocean.  Based  on  the  phenomenology  of  nuclear  detonations  both 

underwater  and  in  the  atmosphere,  the  depth/altitude  axis  is  divided  into  four 

1/3 

regions:  (a)  underwater  (detonation  depth  >  21W  feet,  (b)  atmospheric 

1/3 

(detonation  altitude  >  15W  '  feet),  (c)  near-surface  in  air  (detonation 
1/3 

altitude  <  15W  feet),  and  (d)  near-surface  underwater  (detonation  depth 
<  21W1/3  feet). 

For  the  underwater  case,  the  determination  of  yield  vs  depth  requires  the 
employment  of  one  (or  more)  of  the  extant  Navy-approved  propagation  prediction 
models.  Predictions  are  made  for  the  detonation  s ite-to-recei ver  great  circle 
path  for  assumed  source  placements  throughout  the  water  column.  Yield 
determinations  are  made  based  on  the  predicted  propagation  loss  values  and  the 
known  received  signal. 

Determination  of  yield  vs  altitude  for  the  atmospheric  case  is  based  on  a 
family  of  thirty-six  normalized  relative  received  level  (NRRL)  curves  obtained 
for  the  subject  environment.  This  curve  set  is  obtained  by  the  application  of 
a  simple  scaling  procedure  to  a  corresponding  permanently  stored  set  of  NRRL 
curves  derived  for  a  standard  ocean/ atmosphere  construct.  This  scaling 
involves  knowledge  of  the  maximum  inclination  angle  at  the  ocean  surface  -  in 
the  vicinity  of  the  source  -  associated  with  the  long  range  propagation  modes 

(®s),  and  the  maximum  inclination  angle  at  the  depth  of  minimum  source-to- 
receiver  propagation  loss  (reference  depth)  associated  with  the  long  range 

modes  (9$c)  -  also  in  the  vicinity  of  the  source.  Scaling  is  madp  possible 
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by  the  strong  dependency  exhibited  in  the  equation  for  determination 

of  relative  received  level  {RRL).  The  classification  system  employs  an 

indirect  method  for  determination  of  o  and  o  based  on  the  predicted 

s  sc 

source-to-receiver  propagation  loss  for  both  a  surface  source  (no  upgoing 
energy)  and  a  source  located  at  the  reference  depth.  The  method  makes  the 
assimption  of  cylindrical  spreading  of  the  shot-generated  acoustic  energy  in 
the  vicinity  of  the  receiver.  The  yield-vs-al  titude  curve  for  the  atmospheric 
case  is  defined  by  the  yiel  d/ al  t itude  pairings  within  the  NRRL  curve  family 
for  which  the  value  of  NRRL  is  that  computed  for  the  reference  depth  yield. 

As  for  the  atmospheric  case,  the  method  for  determination  of  yield  vs 
altitude  in  the  near-surface  region  is  based  on  a  set  of  permanently  stored 
NRRL  curves.  In  the  near-surface  region  NRRL  is  assumed  to  be  environmentally 
independent.  The  method  used  for  determining  yield  vs  altitude  using  the 
near-surface  NRRL  curves  is  identical  to  that  employed  for  the  atmospheric 
case . 

In  the  near-surface  underwater  region,  yield  vs  depth  is  determined  by 
simple  linear  interpolation  between  the  yield  computed  for  the  shallowest 
point  on  the  underwater  yield-vs-depth  curve  and  that  computed  for  a  source 
located  right  at  the  surface. 

The  determination  of  NRRL  for  the  reference  depth  is  based  on  the 
passband  filtered,  detected  and  averaged  output  of  an  autocorrelation 
processor  for  the  receiving  hydrophone.  The  required  measurements  are  the 
processor  output  level  due  to  the  received  acoustic  signal  energy  (Sgyj)  and 
the  output  level  due  to  the  antient  noise  energy  just  before  reception  of  the 
shot-gener  ated  signal  (Ngy-j-)*  The  signal  processing  analysis  is  based  on 
the  assumption  that  within  the  operational  frequency  passband  of  the  receiver 
(10  to  30  Hz)  the  spectrin  slope  of  the  detonation  signal  is  -6  dB/octave. 
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This  assumption  is  believed  to  be  valid  for  most  yield/altitude  cont) inations 
for  which  detections  can  be  made.  Based  on  an  averaging  time  of  15  seconds 
and  a  detection  index  of  23  dB,  it  is  estimated  that,  for  the  10  to  30  Hz 
band,  the  minimum  detectable  NRRL  at  a  range  of  2000  nm  in  the  standard  ocean 
environment  is  about  -90  dB.  Due  to  the  large  dynamic  range  requirement  for 
the  system  (~ 120  dB)  an  accurate  automatic  gain  control  (AGC)  must  be  employed 
at  the  receiving  sensor  output.  The  standard  deviation  of  the  error  in  the 
reference-depth  NRRL  determination  is  estimated  to  be  about  3  dB  for  a 
calibrated  linear  receiving  system  and  about  4  dB  for  an  uncalibrated  system. 

For  an  air  detonation  the  determination  of  yield  (or  altitude)  requires 
that  shot  altitude  (or  yield)  be  determined  by  some  other  means.  The 
determination  will  have  a  random  error  that  is  a  nonlinear  function  of  both 
the  error  in  NRRL  and  the  error  in  altitude  (or  yield).  For  an  underwater 
detonation,  a  determination  of  yield  can  be  made  with  a  relatively  high 
accuracy  -  even  if  the  depth  is  not  known  -  because  of  the  small  variation  of 
NRRL  over  the  water  column. 

Other  possible  sources  of  error  in  the  determination  of  yield  vs  depth/ 
altitude  are  those  arising  from  inaccuracies  in  the  derivation  of  the  NRRL 
curves  for  the  standard  ocean/ atmosphere,  and  from  inaccuracies  in  the 
indirect  method  employed  for  determination  of  the  angles  used  in  scaling  NRRL 
to  subject  environments.  These  have  not  been  addressed  in  this  project. 


47 


4.  RECOMMENDATIONS 


The  algorithmic  design  of  a  system  for  the  detection,  localization  and 
classification  of  a  possible  nuclear  detonation  in  the  open  ocean  has  been 
described  in  this  Phase  1  report.  It  is  recommended  that  a  Phase  II  effort  be 
undertaken  for  the  software  implementation  of  the  four  different  yield-versus- 
depth/ altitude  algorithms  (underwater,  atmospheric,  near-surface  in-air,  and 
near-  surface  underwater),  and  their  integration  with  existing  detection  and 
localization  software  into  a  complete  computer  program.  Any  additional 
software  development  needed  by  the  facility  responsible  for  using  this  system, 
such  as  broadband  and  narrowband  signal  processors  and  localization  techniques 
using  time  delays  between  direct  arrivals  and  reverberant  returns  from  such 
ocean  features  as  continental  slopes,  seamounts  and  island  chains  is  also 
recommended  for  Phase  II. 

The  validity  of  the  physical  model  used  to  obtain  the  atmospheric  yield 
and  altitude  functional  dependencies  upon  NRRL  has  not  been  fully  confirmed 
experimentally.  The  derivation  is  based  upon  the  use  of  the  shock  wave  over¬ 
pressure  instead  of  the  ambient  pressure  in  calculating  the  speed  of  sound  and 
density  in  air  at  the  air-water  boundary,  and  the  employment  of  these  values 
in  linear  acoustic  equations.  Qualitative  comparisons  based  upon  nuclear 
detonations  and  observations  from  very  low  yield  chemical  explosive 
experiments  tend  to  support  this  physical  assumption.  We  believe  this  model 
to  be  valid.  If  the  shock  overpressure  assumption  is  correct,  there  will 
exist  a  lateral  wave  whenever  the  overpressure  is  sufficiently  large  that  the 
speed  of  sound  in  air  is  greater  than  the  speed  of  sound  in  water  at  the 
air-water  boundary.  Scale  model  chemical  explosive  experiments  in  air  can  be 
used  to  determine  whether  or  not  a  lateral  wave  exists.  It  is  recommended 
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that  these  tests  be  conducted  in  order  to  validate  the  physical  basis  for  the 
yield-vs-al  t  itude  curves  that  have  been  derived.  These  tests  need  not  impact 
on,  nor  be  a  part  of  the  Phase  II  effort  since  the  yie ld-vs-al titude  curves 
can  be  easily  modified  due  to  any  subsequent  refinements  in  the  physical 
theory. 

Implementation  of  this  system  will  improve  the  overall  capability  of  the 
Federal  Government  to  more  quickly  and  accurately  determine  whether  or  not  an 
at-sea  detonation  has  occurred  and  to  determine  its  location  and  yield-vs- 
depth/ al titude  curve.  This  will  provide  a  greater  opportunity  for  an  aircraft 
to  be  dispatched  to  the  site  and  collect  radioactive  air-borne  particles 
before  they  have  become  too  dispersed  to  be  detected. 

Implementation  will  requi  re  propagation  loss  models,  environmental  data 
banks  and  localization  algorithms,  all  of  which  reside  at  NORDA.  It  is 
anticipated  that  the  integrated  system  will  be  developed  and  implemented  by 
working  closely  with  NORDA  personnel  who  are  familiar  with  the  software 
implementation  of  these  models,  data  banks  and  algorithms.  The  completed 
system  will  reside  at  NORDA  so  that  a  quick  response  can  be  achieved  for  any 
suspected  at-sea  nuclear  detection. 


APPENDIX  A 


DETONATION  LOCALIZATION 

A. 1  Background 

Since  1973,  the  Navy,  under  the  auspices  of  the  Defense  Nuclear  Agency 
(ONA),  has  been  heavily  engaged  in  the  study  of  the  characteristics  of  ocean 
basin  reverberation  generated  by  large  underwater  detonations.  The  overall 
purpose  of  the  effort  is  the  determination  of  the  potential  impact  of 
underwater  nuclear  detonations  on  the  performance  of  underwater  acoustic 
detection  systems.  The  work  has  included  theoretical  studies,  reverberation 
model  development  for  prediction  of  pressure-time  histories,  model  validation 
studies  involving  comparison  of  model  predictions  with  actual  data,  and  the 
planning  and  performance  of  field  experiments  involving  the  detonation  of 
large  chemical  charges  in  the  open  ocean.  Two  computer-implemented 
reverberation  prediction  models  have  been  developed  and  validated  during  this 
program  -  one  by  the  Naval  Surface  Weapons  Center  (NSWC)  and  the  other  by 
Underwater  Systems,  Inc.  (USI),  (as  of  1  March  1984,  Underwater  Systems,  Inc. 
has  become  Underwater  Systems  Group  (USG),  a  division  of  Mandex,  Inc).  The 
major  results  of  the  reverberation  program  are  reported  in  References  4  to  6. 

Additional  investigations  and  analyses  have  been  performed  addressing  the 
question  of  the  feasibility  of  determining  the  location,  yield  and  depth  of  an 
underwater  nuclear  blast  through  analysis  of  the  reverberant  returns  at  a 
single  underwater  sensor,  Reference  7.  The  results  of  this  work  have  been 
quite  encouraging. 

A . 2  Source  Location  Determination 

In  general,  ocean  basin  reverberation  generated  by  underwater  detonations 
can  be  characterized  as  comprising  three  basic  components:  (1)  a  direct  signal 
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of  very  high  amplitude  arriving  via  the  great  circle  track  connecting  the 
source  and  receiver;  (2)  intermittent  large-amplitude  scattered  returns  from 
prominent  basin  reflectors  such  as  coninental  slopes,  island  chains  and 
seamounts;  and,  (3)  lower  level  reverberation  caused  by  scattering  at  the 
ocean  surface  and  basin  bottom  as  well  as  by  multiple  reflections  from  the 
basin  walls.  The  reverberation  pressure-time  history  (levels,  duration  and 
envelope  shape)  at  any  given  site  is  a  strong  function  of  source  level,  source 
and  receiver  locations  within  the  basin,  and  the  basin  geometry.  Figure  A-l, 
taken  from  Reference  6,  shows  an  example  of  a  predicted  100  Hz  pressure-time 
history  for  a  hypothetical  underwater  1-megaton  detonation  located  in  the 
North  Pacific  Ocean.  The  source  and  receiver  locations  are  specified  on  the 
figure.  Also  identified  are  the  basin  reflectors  associated  with  many  of  the 
discrete  returns.  The  ambient  100  Hz  noise  level  was  taken  as  80  dB//yPa/Hz. 
Note  that  for  this  case  reverberant  signals  of  significant  level  are  predicted 
up  to  nearly  two  hours  after  the  reception  of  the  direct  arrival. 

It  is  the  presence  of  the  intermittent  returns  in  large  charge 
reverberation  time  histories  which  afford  the  opportunity  for  estimation  of 
shot  location.  Consider  the  illustration  shown  as  Figure  A-2.  The  figure 
depicts  the  direct  arrival  propagation  track  and  the  propagation  track 
associated  with  the  scattered  returns  from  a  single  prominent  basin  reflector, 
or,  "hot  spot."  The  reverberation  track  comprises  two  legs  -  the  propagation 
path  from  the  source  to  the  reflector  and  the  path  from  the  reflector  to  the 
receiver.  The  ranges  associated  with  the  direct  arrival  track  and  the  first 
and  second  legs  of  the  reverberation  track  are  given  as  Rp,  R^  and  R^ 
respectively.  If  we  assume  some  representative  sound  speed  (c)  for  the  basin, 
the  elapseu  time  after  the  direct  signal  (aT)  at  which  the  reverberation 
return  arrives  at  the  receiver  is  given  by 
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Figure  A-l.  100  Hz  Reverberation  for  a  Hypothetical  1-Meqaton  Detonation  in  the 
North  Pacific  Ocean. 


Figure  A-2.  Geometry  Associated  with  Reverberant  Return  from 
Basin  Reflector  ("Hot  Spot").  (See  text.) 


Assume,  now,  that  a  reverberation  time  history  has  been  obtained  for  an 
underwater  event  and  that  the  location  of  the  detonation  is  unknown.  Also, 
assume  that,  based  on  previous  knowledge  of  the  reverberation  characteristics 
of  the  basin  (obtained,  say,  by  field  experimentation  or  model  predictions) 
that  a  number  of  basin  "hot  spots",  have  been  determined  to  be  associated  with 
specific  peaks  in  the  history.  For  each  "hot  spot",  then,  aT  and  R^<.  R  are 
known.  For  each  case  the  locus  of  points  of  source  locations  which  satisfy 
Eq.  (A-l)  can  now  be  determined.  Neglecting  all  error,  the  resulting  curves 
will  intersect  at  the  true  source  location.  In  actual  practice,  various 
timing  errors  will  result  in  curve  crossings  over  some  finite  area  in  the 
vicinity  of  the  source.  The  nominal  location  is  taken  as  some  centrally 
located  point  within  the  region  of  uncertainty.  The  proper  identification  of 
"hot  spots"  with  specific  reverberant  returns  is,  of  course,  crucial  to  the 
method.  Once  an  estimate  of  a  detonation  site  location  has  been  made,  the 
location  can  be  tested  by  making  a  time-history  prediction  for  the  site  using 
one  of  the  extant  reverberation  prediction  models  and  comparing  the  results 
with  the  measured  history. 

A. 3  Reverberation  Models 

(NOTE:  With  minor  modification,  the 
following  discussion  was  taken  directly 
from  Reference  4. ) 

Two  computer  models  for  predicting  ocean  basin  reverberation  have  been 
developed  and  continually  refined  over  the  past  several  years.  One  model  was 

developed  by  NSWC  and  the  other  was  developed  by  USI  under  contract  to  NSWC. 
Both  models  were  initially  developed  to  predict  reverberation  in  the  North 


Atlantic  Basin.  Experimental  measurements  had  been  made  in  that  area,  and  the 
data  could  be  used  for  validation  purposes.  Both  models  approach  the 
prediction  of  reverberation  in  a  similar  manner,  but  with  some  essential 
differences.  The  NSWC  model  is  designed  more  for  quantity  production,  so 
simplifying  assumptions  have  been  made  to  reduce  computer  running  time.  The 
USI  model  is  basically  a  research  tool  used  to  investigate  special  aspects  of 
reverberation  in  a  basin.  A  brief  discussion  of  each  model  follows. 

A. 3.1  NSWC  Reverberation  Model 

The  Naval  Surface  Weapons  Center  model  is  based  on  the  following 
assunptions  and  procedures: 

a.  Sound  travels  over  long  distances  via  the  sound  channel  along  great 
circle  paths. 

b.  The  basin  boundary  is  defined  by  an  appropriate  depth  contour.  For 
the  North  Atlantic  and  North  Pacific,  a  depth  of  1000  fathoms  has 
been  used . 

c.  Individual  reflectors  in  the  basin,  such  as  islands  or  seamounts,  are 
defined  at  the  same  depth  as  the  basin  boundary. 

d.  Sound  energy  spreads  spherically  to  a  range  of  10  nm  from  the  source 
(or  reflector),  and  cylindrically  thereafter. 

e.  A  modified  Thorp  equation  is  used  to  calculate  the  absorption  in 
water. 

f.  A  constant  scattering  strength  is  assumed  throughout  the  basin.  The 
strength  of  each  reflector  is  determined  by  multiplying  this  constant 
(determined  empirically)  by  a  function  of  slope  angle. 

g.  Individual  signals  are  spread  in  time,  1  second  for  each  100  nm  of 
travel,  to  account  for  multipath  effects.  In  addition,  they  are 
spread  to  account  for  the  size  of  the  reflector. 
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h.  Surface  reverberation  is  provided  by  the  USI  model. 

i.  Multiple  bounces  off  the  basin  walls  may  occur  and  can  be  treated. 

j.  Modeling  a  single  frequency  is  adequate  to  define  the  basin 
reverberation  in  the  10  to  250  Hz  band.  It  has  been  determined 
experimentally  that  only  absorption  corrections  need  be  applied  to 
determine  the  reverberation  at  other  frequencies  within  the  band. 

Making  use  of  these  assumptions  and  procedures,  the  model  calculates  a 
reverberation  record  in  a  straightforward  manner. 

Figure  A- 3  shows  a  simple  basin  model.  The  continental  boundary  and 
obstructions  within  are  defined  by  th  open  circles  and  the  great  circle 
segments  which  connect  them.  Each  segment  that  connects  open  circles  is 
considered  a  reflector.  Boundary  points  which  are  not  connected  define  non¬ 
reflecting  or  open  segnents.  The  source  and  receiver  are  designated  by  the 
points  marked  S  and  R. 

With  the  geometry  defined,  the  computer  model  checks  all  possible  paths 
from  the  source  to  the  receiver.  Of  the  sample  paths  shown,  the  direct  path 
(1)  and  a  single-bounce  path  (5)  reach  the  receiver.  If  desired,  double-  and 
triple-bounce  paths*  are  also  computed.  For  each  path  that  reaches  the 
receiver,  signal  amplitude,  travel  time  and  duration  are  computed. 

A  reverberation  record  is  generated  by  energy  summing  all  the  individual 
arrivals  in  10-second  bins.  These  discrete  arrivals  are  added  to  the  surface 
reverberation  (calculated  frem  the  USI  model)  and  the  anbient  noise. 


*A  double-bounce  path  is  one  that  interacts  with  two  boundary  segments  before 
reaching  the  receiver;  a  triple-bounce  path  interacts  with  three  boundary 
se  gnen  ts . 
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A.  3. 2  US  I  Reverberation  Model 

The  USI  reverberation  model  was  developed  concurrently  with  the  NSWC 
model.  The  model  treats  surface  reverberation  and  boundary  or  seamount 
reflections  with  attenuation  as  a  function  of  grazing  angle.  Multiple 
boundary  bounces,  however,  cannot  be  accommodated. 

The  reverberation  level  due  to  sea  surface  reverberation  is  based  on  the 
following  equation: 

RLs  =  SL  +  SS  +  10  log  A  -  PL1  -  PL2  -  a(Rj  +  R?)  (A-2) 

where 

RLs  is  the  surface  reverberation  level  in  dB//yPa/Hz 
SL  is  the  source  level  in  dB //viP a/H z 
SS  is  the  scattering  strength  of  the  surface  in  dB 
A  is  the  en sonified  surface  area 

PLj  is  the  spreading  loss,  in  dB,  from  the  source  to  the  surface  area 
( 1/R2  to  10  nm,  then  1/R) 

Pl>2  is  the  spreading  loss,  in  dB,  from  the  surface  area  to  the  receiver 
a  is  the  absorption  coefficient 

Figure  A-4  shows  the  geometry  associated  with  the  implementation  of  Eq.  (A-2). 
The  elliptical  ring  depicts  the  area  on  the  surface  from  which  scattered  sound 
arrives  at  the  receiver  if  a  one  second  pulse  width  is  assumed.  The  source 
(S)  and  the  receiver  (R)  are  at  the  foci  of  the  ellipse.  Calculations  are 
made  by  dividing  the  elliptical  ring  into  n  increments  of  angular  width  a©n 
with  each  annular  section  having  an  incremental  area  aAn.  The  scattered 
energy  from  each  aAn  is  then  integrated  over  the  entire  ring. 

Surface  reverberation  versus  time  is  calculated  by  repeating  this  process 
for  each  second  of  time  in  ever-growing  elliptical  rings.  When  the  elliptical 
ring  crosses  a  boundary  segment,  the  angular  increment  subtended  by  the 


boundary  segment  is  excluded  from  the  total  surface  reverberation  level  of 
that  and  all  subsequent  elliptical  rings.  The  calculations  continue  until  the 
ellipse  has  crossed  all  boundary  segments,  or  2  hours  after  the  direct  arrival. 

In  the  coastal  reverberation  part  of  this  model,  a  sloping  continental 
shelf  is  described  by  a  series  of  straight  line  segments  from  which  multiple- 
bounce  reverberation  occurs,  as  illustrated  in  Figure  A-5.  For  each  coastal 
segment,  a  ray  analysis  is  used  to  determine  the  angle  and  depth  of  the 
incident  ray.  Using  appropriate  bottom  loss-  and  bottom  backscattering-ver sus 
grazing-angle  values,  sets  of  reverberation  strengths  and  reverberation  times 
are  calculated  for  the  bottom  encounters  depicted  in  Figure  A-5.  The  times 
are  calculated  relative  to  the  first  bottom  bounce.  When  a  range  of  incident 
angles  and  depths  is  possible,  Idle  reverberation  strengths  and  times  used  are 
an  average  of  the  possible  values. 

The  reverberant  energy  from  each  bottom  encounter  is  distributed  evenly 
over  a  time  interval  appropriate  to  the  size  of  the  boundary  segment.  The 
ranges  from  the  source  to  the  first  encounter  (Rj)  and  from  the  first 

encounter  to  the  receiver  (R^)  are  used  to  determine  spreading  loss  and 
attenuation  for  the  set  of  bottom  encounters.  The  arrbient  noise  level  is 
added  to  the  computed  reverberation  levels. 

The  total  reverberation  record  is  constructed  by  energy  summing  the 
coastal  reverberation  levels  from  all  of  the  boundary  segments  with  the 
surface  reverberation  levels  previously  calculated,  integrated  over  10-second 
interv  al  s. 
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Figure  A-5.  Coastal  Reverberation  Model  -  Upslope  Propagation. 


APPENDIX  8 


SEISMIC  SIGNALS 


The  efficiency  with  which  underwater  detonations  generate  seismic  signals 
detectable  at  teleseismic  ranges  is  well  established  (References  8  and  9). 

One  would  expect  that  atmospheric  detonations  over  the  ocean  would  have  a  very 
low  efficiency  of  seismic  excitation. 

A  first  estimate  of  the  reduction  in  energy  available  for  seismic 
excitation  can  be  obtained  by  considering  the  downward  refraction  of  energy  as 
given  by  the  equations  of  Reference  10.  This  loss  is  given  by: 


L  =  10  log 


(P]c1  +  pc)2  C]2 


4pCp^c^ 


decibels 


(8-1) 


where 


P  *  Pr 


where 


7  +  6p/Pr 

J  +  P/P, 


(8-2) 


p  is  the  density  of  air  behind  the  shock  front  at  the  air-water 
interface 

p0  is  the  ambient  density 

p  is  the  peak  overpressure  (behind  the  shock  front) 

P0  is  the  ambient  pressure  (ahead  of  the  shock) 
p0  *  0.  001253  g/cm3 


PQ  *  14.7  psi 


where 


1/2 


(8-3) 


c  is  the  shock  velocity  in  air  at  the  air-water  interface 


62 


cQ  is  the  ambient  speed  of  sound  (ahead  of  the  shock  front) 
p  is  the  peak  overpressure  (behind  the  shock  front) 

PQ  is  the  ambient  pressure  (ahead  of  the  shock) 
cQ  =  33528cm/ sec  (11 00ft/ sec) 

PQ  =  14.7  psi 

pj  is  the  density  of  water  behind  the  shock  front  at  the  air-water 
interface 

Cj  is  the  shock  velocity  in  water  at  the  air-water  interface 

These  equations  are  consistent  with  what  we  have  done  in  analyzing 
the  acoustic  problem.  The  shock  wave  velocity  and  density  in  air  is  used 
in  linear  acoustic  equations  which  properly  handle  refraction  and 
sphericity  to  predict  the  downward  propagating  signal.  The  results  are 
shown  in  Figure  B-l  which  shows  the  loss  as  a  function  of  altitude  for 
parametric  yields  ranging  from  0.01  kiloton  to  10  megatons.  Note  that 
there  is  a  yield-altitude  region  for  which  the  value  of  L  is  negative. 
That  is,  there  is  a  gain  rather  than  a  loss.  This  is  caused  by  the 
downward  refraction  of  the  energy  which  is  sufficient  to  overcome  the 
impedance  mismatch  at  high  enouc^i  shock  velocities. 

The  loss  factor  L  describes  only  the  loss  of  energy  transmitted 
throu^i  the  interface.  As  compared  to  an  underwater  detonation,  there 
are  additional  losses.  This  includes  3  dB  corresponding  to  the  loss  of 
the  upward  propagating  energy  (reflected  from  the  sea  surface)  which 
exists  for  the  underwater  detonation,  and  about  6  dB  corresponding  to  the 
loss  of  the  bubble  pulse  energy  which  is  a  major  contributor  to  the  low 
frequencies  considered  in  seismic  signals.  Taking  these  factors  into 
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Loss  (dB) 

Energy  Loss  Through  the  Air-Water  Interface  as  a 
Function  of  Altitude  for  the  Indicated  Yields. 


Figure  B-l 


account,  the  results  of  Figure  B-l  can  be  replotted  as  shown  in  Figure 
B-2. 

The  solid  lines  of  Figure  B-2  show  the  yie  Id/ al  titude  curves  for 
fixed  total  excitation  losses  (parameterized).  The  dashed  lines  provide 
a  means  for  relating  deep  underwater  and  atmospheric  detonations.  For 
example,  a  200-kiloton  atmospheric  detonation  at  an  altitude  of  about 
240  feet  will  generate  the  same  seismic  signal  as  a  deep  underwater 
detonation  of  250  kilotons.  Similarly,  a  20-kiloton  detonation  at  an 
altitude  of  110  feet  is  equivalent  to  a  deep  underwater  detonation  of 
about  25  ki  lotons. 

The  numerical  values  provided  by  Figure  B-2  should  not  be  taken  too 
literally.  Energy  of  atmospheric  detonations  is  crudely  modeled  by 
10  log  W  decibels,  where  W  is  the  yield,  and  the  upward  going  energy  and 
bubble  pulse  of  the  underwater  detonation  are  only  crudely  treated. 
Further,  on  fundamental  considerations,  we  would  expect  that  for  a 
surface  detonation  the  available  energy  would  be  reduced  by  9  dB.  The 
computations  indicate  that  the  efficiency  of  excitation  is  higher  for  the 
low  altitude  atmospheric  detonations. 

However,  it  is  safe  to  conclude  that  low  altitude  nuclear 
detonations  can  generate  detectable  seismic  signals  comparable  to  those 
produced  by  underwater  detonations.  The  parametric  line  designated  by 
-1  dB  (gain)  roughly  represents  the  yield/ altitude  curve  above  which  the 
energy  available  for  generating  seismic  signals  reduces  rapidly. 
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Figure  B-2.  Yield/Altitude  Curves  (Solid  I 
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c  is  the  shock  velocity  in  air  at  the  air-water  interface 
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